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The extended Hubbard model in the zero-bandwidth limit is studied. The effective Hamiltonian 
consists of (i) on-site U interaction and intersite (ii) density-density interaction W and (iii) Ising- 
like magnetic exchange interaction J (between the nearest-neighbors). We present rigorous (and 
analytical) results obtained within the transfer-matrix method for ID-chain in two particular cases: 

(a) W = 0 and n = 1; (b) U -l-oo and n = 1/2 {W yf 0, J 7 ^ 0). We obtain the exact formulas 
for the partition functions which enables to calculate thermodynamic properties such as entropy, 
specific heat (c), and double occupancy per site. In both cases the system exhibits an interesting 
temperature dependence of c involving a characteristic two-peak structure. There are no phase 
transitions at finite temperatures and the only transitions occur in the ground state. 
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I. INTRODUCTION 

The interplay between density-density and magnetic 
interactions is relevant to a broad range of important ma¬ 
terials such as manganites, multiferroics, organics, and 
other strongly correlated electron systems HHg. 

In this paper we present some exact results obtained 
within transfer matrix method for the zero-bandwidth 
extended Hubbard model with density-density and Ising- 
like magnetic interactions on the one dimensional chain 
(d = 1). The ID-Hamiltonian considered has a form 

L 

H = '^ [Ufii^fii^ + Wfi.fii+i - - fifii], ( 1 ) 

where denotes the creation operator of an electron 
with spin a {a =t,4.) at site i, hia = c^^Cia, fn = 
and sf = (1/2) (hi^ — n^). i J- 1 is the nearest neighbor 
of the Tsite in the chosen direction (from two possible 
directions in a chain). We assume the periodic boundary 
conditions, i.e. u-i+io- = rrio-, where L is a number of 
sites in the chain and tiia = {fiia). Jq = zJ, where z = 2 
is a number of the nearest neighbors. 

All the terms of Hamiltonian Q commute with one 
another and are diagonal in the representation of occu¬ 
pancy numbers. It is convenient to use the transfer ma¬ 
trix method HO] to find the grand partition function Z. 
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Hamiltonian 0 can be treated as a simple effective 
model of insulators, in which interactions U, W and J 
are assumed to include all the possible contributions and 
renormalizations. Notice that ferromagnetic (J > 0) in¬ 
teractions are simply mapped onto the antiferromagnetic 
ones (J < 0) by redefining the spin direction on one sub¬ 
lattice in lattices decomposed into two interpenetrating 
sublattices. Thus, we restrict ourselves to a case of J > 0. 

Exact solutions of model (Q for some particular cases 
have been obtained for the one-dimensional case (T > 0) 
employing the method based on the equations of motion 
and Green’s function formalism or the transfer- 

matrix method [HHn]. Extensive mean-field studies (ex¬ 
act result in d —-l-oo) [THl - 127] and some Monte Carlo 
simulations {d = 2) |28f[5D] of model 0 have been also 
performed. Moreover, the exact ground state (T = 0) 
results have been found for 2 < d < -too ISIMl- 

We present rigorous results for partition functions 
obtained within the transfer-matrix method for one¬ 
dimensional model Q in two particular cases: (a) W = 0 
and n = 1; (b) [/ —>■ -too and n = 1/2 {W 7^ 0, J yf 0). 

II. RESULTS AND DISCUSSION 

(a) For the case of W = 0, a typical element of the 
transfer matrix for model Q is defined as 

Pi,i+i = {ni^nix.\P\ni+i^ni+i^) = ( 2 ) 

= exp {-(3 [{U/2) (nit^4 + 

- (/x/2)ni+i - J{ni^ - nil)(n^+l■^ - rn+n)]} , 
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FIG. 1. The specific heat c as a function of ksT/Jo for several values of U/Jo (as labeled); n = 1, W = G. 


where \ni^nii) e {| 00 ), | 01 ), | 10 ), | 11 )} denotes a single¬ 
site state at site i, (3 = \/{kBT) is the inverse tempera¬ 
ture and ks is the Boltzmann constant. One obtains 16 
matrix elements and the problem is reduced to diagonal- 
ization of the matrix P of the form 

/ I Xo Xo uqxI \ 

p _ Xo moxl rriQ^xl uoxj^ , . 

Xo rriQ xg moxg Uoxg 

\ Moa^o ^oa^o uqXq ugx^ / 

where xo = exp(/3/i/2), uq = exp(—/3U/2) and 
mo = exp(/?Jo/2). One can show that three eigen¬ 
values of P (A;, Z = 1, 2, 3) are roots of a cubic equation: 

A^—A^(l -I- 2mox -I- ux^) — A {xmQ ^(1 — mo)x 
X [x -I- mox -h mgX -I-mQ(2-I- X -I-2ux^)]}-I- (4) 

-I- x^mg ^(1 - mo)^(l -h mo)(l -h ux^) = 0, 

where x = Xg and u = Uq, while A 4 = 0 . 

So far the number N of particles in the chain has not 
been specified. It can be done in a standard way by solv¬ 
ing the following equation: N = — (90/9/i)^, where 0 
is the grand canonical potential, 0 = — fc^TlnZ. In the 
thermodynamic limit L —>■ 00 the grand sum of states Z is 
derived as Z = where Am is the maximum eigenvalue 
of P (assumed to be nondegenerate). Therefore, the 
equation for N can be rewritten as OXm/ dx = uXm jx, 
where n = N/L is electron density in a system. 

In the case of half-filling (n = 1), the condition for N 
can be solved analytically for arbitrary U and in such a 
case the chemical potential is derived as ^ = Z7/2 and Am 
takes the form 


where JC — -t- Zi — 4.^2 — 4^3 - 1 - 16.^4 -t- 4.^5 -t- 2 .^g, 
Y = U-Jo, Z^= exp (2/3Jo), Z 2 = exp {-PY/2), 
Z 3 = exp (/3A/2), Z 4 = exp (/3B/2), Z 5 = exp (-PY), 
Zq = exp {j3Jo), and A = 3Jo — U, B = 2Jq — U. 

(b) The limit U —)■ -l-oo corresponds to the subspace 
where the double occupancy of sites is excluded (by elec¬ 
trons for n < 1 or holes for n > 1). For this case the 
transfer matrix elements for model Q are defined as 

Ri^i+i = {ni^nil\R\n^+l^ni+ll) = ( 6 ) 

= exp {-/3 [W{ni-^ + nii){ni+i^ + rn+ii) 

- (m/2) (n, + Tii+i) - J{ni^ - n4)(ni+i^ - n,+i;)]} , 

where G {| 00 ), | 01 ), | 10 )} denotes a single-site 

state at site i in the limit U —>■ -l-cx). Therefore, in this 
case the matrix R has the following form 

( 1 Xo Xo \ 

Xo mox^wo mg ^xgwo , (7) 

Xo mg ^x^wo moxlwo J 

where xq = exp(/3^/2), wg = exp(—/3Wo/2) and 
mo = exp(/3Jo/2). The eigenvalues of P are roots 
A; (I = 1 , 2 ,3) of the following cubic equation: 

A^ — A^(l -t- 2i71oXWo) 

— Ax[2 — 2mo'Wo A xw^ttiq'^ — m^)] A ( 8 ) 

-I- (1 — mg ^)x^w(l -I- mg — 2moWQ^) = 0, 
where x = Xg and w = Wg- 

In this case the equation for N can be 
solved analytically for n = 1/2. One finds that 
/r = ITo/2 — fesTln [2cosh (/IJg/2)], and Am is derived 

Am = 1 + exp {pWg/A) v'sech (,5Jo/2). (9) 


Am = 1 + exp 




+ -exp 



The knowledge of explicit form of the sum of states Z 
allows us to obtain thermodynamic characteristics of the 






















3 



k^T/J^ W/J 


FIG. 2. The entropy s/fc_B = sasafunction (a) offcsT/Jo and (b) oiWj J {U —>■ +oo, n — 1/2, values of other model parameters 
as labeled). 


system for arbitrary temperature. Local magnetic mo¬ 
ment 7 is defined by: 7 = (1/2L) (|ni^ — n^l). It is 

related with the double occupancy D per site (defined by 
the formula: D = (1/L) = {df/dU)rp) by the re¬ 

lation: 'J = n/2 — D, where f = uj + nfj, is the free energy 
of the system per site (w = fl /L). The entropy s and the 
specific heat c (per site) can be derived as: s = —df/dT 
and c = —T{d^f/dT^). Because the explicit forms of the 
partition function Z in both cases are known and the 
derivation of the above thermodynamical characteristics 
(i.e. D, s, c) is rather straightforward (w = —fcsTlnAM), 
below we only summarize the most important conclusions 
following from the analysis of Eqs. ^ and @. 

One can observe that the system exhibits an interesting 
temperature dependence of c involving a characteristic 
two-peak structure for some values of model parameters 
in the cases analyzed. In both cases considered above 
there are no phase transitions at finite temperatures (in 
the agreement with Mermin-Wagner theorem |35| ) and 
the only transitions can occur in the ground state. 

A. The case of VF = 0 and n = 1. 

For large U/ Jq c exhibit two peak structure, whereas 
for U ^ 6 peaks merge and there is a single peak in T- 
dependence of c (labeled as Ti -|-r 2 , cf. Fig. [^. For 
— 1 < U/Jq < —0.9 two peaks of c appear again. The 
broad one (at higher temperature Ti) is connected with 
continuous changes in short-range charge on-site order¬ 
ing (associated with U term). The narrow one (at 
lower temperature T 2 < Ti) is connected with short- 
range intersite magnetic ordering (J term). With de¬ 
creasing of U/ Jq their locations move towards lower 
temperatures. If U/Jq < —1 the single maximum of 


c (connected with short-range on-site ordering) exists 
only and it moves toward higher temperatures with 
increasing of \U\/Jo. In the limit U —>■-boo (n = 1) 
the specific heat exhibits a single peak described by 
= fcs [(/3Jo/2) X sech (/3Jo/2)]^, which corresponds 
to ID-Ising model in the absence of magnetic field (the 
peak connected with on-site ordering is „located” at 
Ti —>■ -boo). The divergence of c at U/Jq —>■ —1 and 
T 2 ^ 0 indicates that the first-order transition occurs 
between the nonordered state of double occupied sites 
(7 = 0) and the ferromagnetic homogeneous phase (sta¬ 
ble for U/Jq>—1), where all sites are singly occu¬ 
pied (7 = 1 / 2 ). One can derive the same conclusion 
of r = 0 properties of the system from a behavior of 
the entropy s. For U/Jq > — 1 the system is magneti¬ 
cally ordered with s(0) = 0 (a number of states is g = 2, 
s = kss = {kB/L)lng). For U/Jq < —1 the system con¬ 
sists of nonordered on-site electron pairs and s(0) = ln2 
{g = 2^). At C//Jo —)■ —1 s(0) = In2. Notice that in the 
limit T —> 00 the entropy s —> 2 In 2 for any U/ Jq. 

Our results for W = 0 and n = 1 are in an agreement 
with the results of Ref. m obtained using the Green’s 
function formalism, whereas the numerical analyses of Q 
and the condition for N (for arbitrary n or fi) should be 
consistent with the results of Ref. m- 

B. The case of n = 1/2 for U -boo. 

At T = 0 for IT/J = 1 the transition between the ho¬ 
mogeneous charge-ordered (CO) phase (for W/J > 1) 
and phase separated (PS) state occurs, cf. also Ref. |24| . 
For W/J<1 the system is divided into two equal-sized 
domains: one ferromagnetically ordered completely filled 
by electrons (n = 1) and the other empty (n = 0). The 
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behavior of c at T > 0 is very similar to that discussed 
in the previous case (simplifying, for qualitative discus¬ 
sion only U/Jq <->■ —WjJ replacement is needed, short- 
range charge order peak in c is associated with W term). 
In the limit W/ J —>■ —(X) c exhibits a single maximum 
described by the characteristic dependence for ID-Ising 
model, but in such a case the specific heat c* of the sys¬ 
tem is 2 times smaller than the result obtained in 
a case of n = 1, W = 0 {c* = j2). It can be derived 

that at r = 0: (i) for W/J > 1: s(0) = (1/2) ln2 (the 
CO phase, g = 2^'^/^+^)) and (ii) for W/J < 1: s(0) = 0 
(the PS:F/NO state, g = 2L). liW = J the ground state 
is highly degenerated and s(0) = ln(l -|- v^) (cf. Fig. [^. 
In the limit T —?> oo the entropy s —>■ (3/2) In 2. 

The detailed discussion of thermodynamic properties 
of one-dimensional model 0 in a general case will be the 


subject of a subsequent paper. 


ACKNOWLEDGMENTS 

The work (K.J.K, S.R.) has been financed by National 
Science Centre (NCN, Poland) as a research project 
under grant No. DEC-2011/01/N/ST3/00413 and 
a doctoral scholarship No. DEC-2013/08/T/ST3/00012. 
K.J.K. thanks for the financial support from the ESF - 
OP “Human Capital” - POKL.04.01.01-00-133/09-00 - 
Proinnowacyjne ksztaicenie, kompetentna kadra, absol- 
wenci przysziosci^\ K.J.K. and W.K. thank the Foun¬ 
dation of Adam Mickiewicz University in Poznan for the 
support from its scholarship programme. 


[1] R. Micnas, J. Ranninger, S. Robaszkiewicz, Rev. Mod. 
Phys. 62, 113 (1990). DOT: 10.1103/RevModPhys.62.113 

[2] T. Goto, B. Liithi, Adv. Phys. 52, 67 (2003). DOT: 
10.1080/0001873021000057114 

[3] E. Dagotto, T. Hotta, A. Moreo, Phys. Reports 344, 1 
(2001). DOT: 10.1016/80370-1573(00)00121-6 

[4] J. van den Brink, D.I. Khomskii, J. Phys.: Con- 
dens. Matter 20, 434217 (2008). DOI: 10.1088/0953- 
8984/20/43/434217 

[51 W.R. Czart, S. Robaszkiewicz, Acta Phys. Pol. A 106, 
709 (2004); 

W.R. Czart, S. Robaszkiewicz, Phys. Status Solidi (b) 
243, 151 (2006). DOI: 10.1002/pssb.200562502 

[ 6 ] W.R. Czart, S. Robaszkiewicz, Acta Phys. Pol. A 109, 
577 (2006); 

W.R. Czart, S. Robaszkiewicz, Material Science - Poland 
25, 485 (2007). 

[7] K. Kapcia, Acta Phys. Pol. A 121, 733 (2012). 

[ 8 ] K. Kapcia, J. Supercond. Nov. Magn. 26, 2647 (2013). 
DOI: 10.1007/sl0948-013-2152-l 

[9] K.J. Kapcia, Acta Phys. Pol. A 127, 204 (2015); DOI: 
10.12693/APhysPolA. 127.204 

K.J. Kapcia, J. Supercond. Nov. Magn. 28, 1289 (2015). 
DOI: 10.1007/sl0948-014-2906-4 

[10] G.F. Newell, E.W. Montroll, Rev. Mod. Phys. 25, 353 
(1953). DOI: 10.1103/RevModPhys.25.353 

[11] F. Mancini, F.P. Mancini, Phys. Rev. E 77, 061120 
(2008). DOI: 10.2478/sll534-012-0017-z 

[12] F. Mancini, E. Plekhanov, G. Sica, Cent. Eur. J. Phys. 
10, 609 (2012). DOI: 10.2478/sll534-012-0017-z 

[13] F. Mancini, E. Plekhanov, G. Sica, Eur. Phys. J. B 86, 
224 (2013). DOI: 10.1140/epjb/e2013-40046-y 

[14] P. Mancini, E. Plekhanov, G. Sica, J. Phys.: 
Conf. Series 391, 012148 (2012). DOI: 10.1088/1742- 
6596/391/1/012148 

[15] R.A. Bari, Phys. Rev. B, 2662 (1971). DOI: 

10.1103/PhysRevB.3.2662 

[16] G. Beni, P. Pincus, Phys. Rev. B 9, 2963 (1974). DOI: 
10.1103/PhysRevB.9.2963 

[17] R.S. Tu, T.A. Kaplan, Phys. Status Solidi (b) 63, 659 
(1974). DOI: 10.1002/pssb.2220630229 

[18] R. Micnas, S. Robaszkiewicz, K.A. Chao, Phys. Rev. B 
29, 2784 (1984). DOI: 10.1103/PhysRevB.29.2784 


[19] K. Kapcia, W. Klobus, S. Robaszkiewicz, Acta. Phys. 
Pol. A 118, 350 (2010). 

[20] K. Kapcia, S. Robaszkiewicz, J. Phys.: Con- 

dens. Matter 23, 105601 (2011). DOI: 10.1088/0953- 
8984/23/10/105601 

[21] K. Kapcia, S. Robaszkiewicz, J. Phys.: Con- 

dens. Matter 23, 249802 (2011). DOI: 10.1088/0953- 
8984/23/24/249802 

[22] K. Kapcia, S. Robaszkiewicz, Acta Phys. Pol. A 121, 
1029 (2012). 

[23] W. Klobus, K. Kapcia, S. Robaszkiewicz, Acta. Phys. 
Pol. A 118, 353 (2010). 

[24] K. Kapcia, W. Klobus, S. Robaszkiewicz, Acta Phys. 
Pol. A 121, 1032 (2012). 

[25] S. Robaszkiewicz, Phys. Status Solidi (b) 59, K63 (1973). 
DOI: 10.1002/pssb.2220590155 

[26] S. Robaszkiewicz, Phys. Status Solidi (b) 70, K51 (1975). 
DOI: 10.1002/pssb.2220700156 

[27] S. Robaszkiewicz, Acta Phys. Pol. A 55, 453 (1979). 

[28] S. Murawski, K. Kapcia, G. Pawlowski, S. 

Robaszkiewicz, Acta Phys. Pol. A 121, 1035 (2012). 

[29] S. Murawski, K.J. Kapcia, G. Pawlowski, S. 
Robaszkiewicz, Acta Phys. Pol. A 126, A-110 (2014). 
DOI: 10.12693/APhysPolA.126.A-110 

[30] S. Murawski, K.J. Kapcia, G. Pawlowski, S. 
Robaszkiewicz, Acta Phys. Pol. A 127, 281 (2015). 
DOI: 10.12693/APhysPolA. 127.281 

[31] S.A. Pirogov, Ya.G. Sinai, Theor. Math. Phys. 25, 1185 
(1975); DOI: 10.1007/BP01040127 

S.A. Pirogov, Ya.G. Sinai, Theor. Math. Phys. 26, 39 
(1976). DOI: 10.1007/BP01420583 

[32] J. Jedrzejewski, Z. Phys. B 48, 219 (1982); DOI: 
10.1007/BP01420583 

J. Jedrzejewski, Z. Phys. B 59, 325 (1985). DOI: 
10.1007/BP01307438 

[33] U. Brandt, J. Stolze, Z. Phys. B 62, 433 (1986). DOI: 
10.1007/BP01303574 

[34] J. Jedrzejewski, Physica A 205, 702 (1994). DOI: 
10.1016/0378-4371(94)90231-3 

[35] N.D. Mermin, H. Wagner, Phys. Rev. Lett. 17, 1133 
(1966). DOI: 10.1103/PhysRevLett.l7.1133 








